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Abstract. We present new high-resolution observations of the nucleus of the counter-rotating LINER NGC 4826, 
5^ ' made in the J=l — and 3=1 — 1 lines of 12 CO with the IRAM Plateau de Bure mm-interferometer(PdBI).The 

CO maps, which achieve 0.8"(16pc) resolution in the 2-1 line, fully resolve an inner molecular gas disk which is 
truncated at an outer radius of 700 pc. The total molecular gas mass (3.1x10 8 Mq) is distributed in a lopsided 
nuclear disk of 40 pc radius, containing 15% of the total gas mass, and two one-arm spirals, which develop at 
different radii in the disk. The distribution and kinematics of molecular gas in the inner 1 kpc of NGC 4826 
show the prevalence of different types of m=l perturbations in the gas. Although dominated by rotation, the 
gas kinematics are perturbed by streaming motions related to the m=l instabilities. The non-circular motions 
associated with the inner m—l perturbations (lopsided instability and inner one-arm spiral) agree qualitatively 
with the pattern expected for a trailing wave developed outside corotation ('fast' wave). In contrast, the streaming 
motions in the outer m=l spiral are better explained by a 'slow' wave. A paradoxical consequence is that the 
inner m=l perturbations would not favour AGN feeding. An independent confirmation that the AGN is not 
being generously fueled at present is found in the low values of the gravitational torques exerted by the stellar 
potential for R<530 pc. The distribution of star formation in the disk of NGC 4826 is also strongly asymmetrical. 
The observed asymmetries, revealed by HST images of the inner disk, follow the scales of the various m=l 
perturbations identified in the molecular gas disk. Massive star formation is still vigorous, fed by the significant 
molecular gas reservoir at R<700 pc. There is supporting evidence for a recent large mass inflow episode in 
NGC 4826. The onset of m=l instabilities of the type observed in NGC 4826 may be a consequence of secular 
evolution of disks with high gas mass contents. 

These observations have been made in the context of the NUclei of GAIaxies (NUGA) project, aimed at the study 
of the different mechanisms for gas fueling of Active Galactic Nuclei (AGN). 
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1. Introduction 

1.1. Feeding Active Galactic Nuclei: the Nuclei of 
Galaxies (NUGA) project 

Although there is observational evidence that the major- 
ity of galaxies contain super-massive black holes in their 
nuclei, the existence of nuclear activity is far from univer- 
sal (Ho et al. 1997). It is commonly accepted that nuclear 
activity results from the feeding of a massive black hole 
by the infall of gas from its host galaxy. However, there 
is no consensus on which mechanisms are responsible for 
removing the angular momentum from the gas and driv- 
ing infall down to scales of tens of pc. Moreover, it is 
unknown whether these mechanisms are at work only in 
active galaxies, or, alternatively, if the key difference be- 
tween active and quiescent objects is the availability of 
gas supply to the nucleus. On large scales, several dynam- 
ical perturbations induced by galaxy collisions, mergers 
and mass accretion (Hcckman ct al. 1986) can effectively 
drive infall of gas to scales of ~ 1 — 3 kpc. Numerical sim- 
ulations have also shown that gravity torques from barred 
potentials are efficient at funneling gas into the center of a 
galaxy. Observational support for the role of bars in driv- 
ing gas inflow at these scales has been found (Sakamoto 
et al. 1999). However, when an Inner Lindblad Resonance 
(ILR) is present, the gas accumulates in rings (Combes 
1988; Piner et al. 1995; Buta & Combes 1996); subse- 
quently there is virtually no gas infall to the nucleus, as 
gravity torques are positive inside the ILR. 

To explain the mass transfer from the ~ 1 kpc scale 
to radii of tens of pc, we must deal with the non-trivial 
problem of angular momentum removal. Possible solutions 
have invoked the onset of non-axisymmetric dynamical 
perturbations such as nested nuclear bars (Shlosman et 
al. 1989; Fricdli & Martinet 1993); lopsidedness or m=l 
instabilities (Shu ct al. 1990; Kormendy & Bender 1999; 
Garcfa-Burillo et al. 2000); warped nuclear disks (Pringlc 
1996; Schinnerer et al. 2000a, 2000b); and nuclear spiral 
density waves (Englmaicr & Shlosman 2000). 

The study of interstellar gas in the nuclei of galax- 
ies is a fundamental tool for understanding nuclear ac- 
tivity and its relation to circumnuclear star formation. 
From the theoretical point of view, there is an increasing 
body of evidence that purely gaseous density waves (spi- 
rals, bars, warps or lopsided instabilities) may be driving 
gas infall to the Active Galactic Nucleus (AGN) (Heller & 
Shlosman 1994; Elmcgrcen ct al. 1998; Regan & Mulchaey 
1999). Within the central kiloparsec, most of the gas is in 
the molecular phase, which makes CO lines the best trac- 
ers of nuclear gas dynamics. Up to now, CO surveys of 
galaxies made with single-dish telescopes were hampered 
by insufficient spatial resolution (Kcnncy & Young 1988; 
Hcckman et al. 1989; Young et al. 1995; Braine et al. 1993; 
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Casoli et al. 1996; Vila-Vilaro et al. 1998). Most CO inter- 
ferometer surveys of nearby spirals (Sakamoto et al. 1999; 
Regan et al. 2001; Heifer et al. 2003) mapped CO(1-0) disk 
emission at low spatial resolution (4"-7") and moderate 
sensitivity (dctcctability thresholds T, gas >180Mq/pc 2 ). 
Furthermore, until very recently (Jogee et al. 2001), the 
published survey samples have included very few AGN. 

A study of the gas fueling in AGN absolutely requires 
high (~0.5"-l") spatial resolution, given the small (~50- 
100 pc) linear scales involved. We also need high sensitivity 
to provide a high dynamic range in the synthesized maps. 
The NUclei of GAlaxies project-NUGA-(Garci'a-Burilk) 
et al. 2003) is a high-resolution, high-sensitivity CO survey 
of a sample of 12 nearby AGN which spans the sequence 
of activity types (Seyfertl, Seyfert2 and LINERs). The 
survey is being carried out with the IRAM Plateau de 
Bure mm-interferometer (PdBI) in France, which offers 
the best combination of sensitivity and resolution crucial 
for this project. NUGA aims at reaching spatial resolu- 
tions <1" and will allow the acquisition of CO maps with 
high dynamic ranges for most of the targets. Our objective 
is to determine the distribution and dynamics of molec- 
ular gas in the inner 1 kpc of the nuclei with resolutions 
of ~10-50pc, and to study systematically the different 
mechanisms for gas fueling of the AGN. 

The NUGA project relies on a multi-wavelength ap- 
proach: the sample of galaxies observed at PdBI has 
been defined based on the availability of high-quality 
optical and near-infrared (NIR) images. These images 
have been obtained with both ground-based telescopes 
and the Hubble Space Telescope (HST). The optical and 
NIR counterpart of NUGA will provide the stellar poten- 
tials, the star formation history and the dust distributions 
for the galaxies in our sample. We will thus be able to 
study quantitatively the debated AGN-starburst connec- 
tion. NIR maps will be used to compute the stellar poten- 
tials. These will be the basis for self-consistent numerical 
simulations of the gas dynamics (obtained from CO) in 
real case scenarios (defined by the optical and NIR im- 
ages). The long term aim is to complete a super-sample 
of 25-30 objects observed by consortium members, within 
and outside NUGA, with the IRAM array. The initial pa- 
pers will exploit on a case-by-case basis the data obtained 
in the context of the NUGA project. 

1.2. The counter-rotator LINER NGC4826 

In this first paper, we study the distribution and dy- 
namics of molecular gas in the nucleus of the LINER 
NGC4826, using high-resolution (0.8"-3") observations 
made in the 1-0 and 2-1 lines of 12 CO with the IRAM 
array. NGC4826, also known as the 'Black Eye' or 'Evil 
Eye galaxy' due to its optical appearance, hosts two 
nested counter-rotating gas disks of comparable mass. 
The inner disk extends to a radius of ^50" and contains 
-10 7 M Q of atomic (HI) gas (Braun et al. 1992, 1994) 
and — 2.3x10 8 Mq of molecular (H 2 ) gas (Casoli & Gerin 
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Fig. 1. 12 CO(1-0) velocity-channel maps observed with the PdBI in the nucleus of NGC 4826 with a spatial resolution 
of 3"x2.4" at PA=44° (beam is plotted as a filled ellipse in the bottom left corner of each panel). We show a field 
of view of 75", i.e. ~1.8 times the diameter of the primary beam at 115 GHz. The phase tracking center is indicated 
by a cross at aj2ooo=12 /l 56 m 43.88 s and <5j2ooo=21°41'00.1". Velocity-channels are displayed from v=-190 km s _1 to 
v=190 km s -1 in steps of 10 km s _1 . Velocities are in LSR scale and referred to v=v o =408 km s -1 . Contour levels are 
-3<t, 3a, 10a, 20cr, 40a and 78a where the 1-sigma rms a=4.0 mJy beam -1 . There are very few pixels with negative 
fluxes below -3a inside the field-of-view. 



1993). The outer disk, extending from ~80" to ~9.8', ro- 
tates in the opposite sense to the inner gas and contains 
-1O 8 M of HI (Braun et al. 1992, 1994). Rix et al. (1995), 
by studying the stellar kinematics along the principal axes 
of NGC 4826, found that the stars rotate at all radii with 
the same sense as the inner disk providing strong evidence 
that stars and gas are coplanar. However, the kinematics 
of the ionized gas, analyzed by Rubin (1994) and Rix et 
al. (1995), are considerably more complex. Within radii 
<30", ionized gas co-rotates with HI and the stars, fol- 
lowed by a strong kinematic disturbance in the region 
30" <r<100", where rotation velocity is close to zero. 
Most remarkably, Rix et al. (1995) found evidence for ra- 
dial inflow along the minor axis. 

Casoli & Gerin (1993) made low-resolution 12 CO(l- 
0) (22") and 12 CO(2-l) (12") maps of NGC4826's disk 
using the IRAM 30m telescope. Their maps showed that 
the co-rotating molecular disk ends quite abruptly at a 



radius of 45-50", i.e., at a significantly small extent rel- 
ative to the optical galaxy size (1)25=10'). These low- 
resolution maps could not reveal the small-scale structure 
and kinematics of the compact molecular disk. NGC 4826 
was also observed as a part of the interferometer surveys of 
OVRO (Sakamoto et al. 1999) and BIMA-SONG (Heifer 
et al 2003). These interferometer maps confirmed that the 
molecular gas disk of NGC 4826 extends out to a radius 
of 700 pc. The OVRO map of NGC 4826, published by 
Sakamoto et al. (1999), shows tantalizing evidence of an 
asymmetrical distribution of molecular gas in the inner 
1 kpc. 

The maps presented in this paper fully resolve the 
molecular gas distribution in NGC 4826 out to a radius 
of r=35"(700pc). The high spatial resolution (0.8"=16pc 
in the 2-1 line) and high dynamic range (~70-100) of 
these observations allow the study of the complex pat- 
tern of gravitational instabilities at work from lkpc to 



4 



S. Garci'a-Burillo et al.: Molecular Gas in NUclei of GAlaxies (NUGA) 




Fig. 2. Same as Fig.Q]but for the 2-1 line of 12 C0. Spatial resolution reaches l.l"x0.8" at PA=34° (beam is plotted 
as a filled ellipse in the bottom left corner of each panel). We show a field of view of 43", i.e. ~2 times the diameter 
of the the primary beam at 230 GHz. Velocity-channels are displayed from v=-172.5 km s _1 to v=167.5 km s _1 in 
steps of 10 km s _1 , with same reference as used in Fig. ^ Contour levels are -4a, 4a, 8a, 16a, 32a and 62a, where 
the 1-sigma rms a— 5.7 mJy beam" 1 . There are very few pixels with negative fluxes below -4a inside the field-of-view. 

urations of the array (Guilloteau et al. 1992). We ob- 
served simultaneously the J=l-0 and J=2-l lines of 12 CO 
in a single field centered at aj 2 ooo=12' l 56 m 43.88 s and 
<5/20oo=21°41'00.1"; the primary beam size is 42" (21") in 
the 1-0 (2-1) line. The spectral correlator was split in two 
halves centered at 115.114 GHz and 230.224 GHz, respec- 
tively, i.e., the transition rest frequencies corrected for an 
assumed recession velocity of v o (iS'i?)=408km s _1 . The 
correlator configuration covers a bandwidth of 580 MHz 
for each line, using four 160MHz-wide units; this is equiv- 
alent to 1510km s~ 1 (755km s^ 1 ) at 115GHz (230GHz). 
The correlator was regularly calibrated by a noise source 
inserted in the IF system. Visibilities were obtained using 
on-source integration times of 20 minutes framed by short 
(~ 2 min) phase and amplitude calibrations on the nearby 
quasars 1308+326 and 3C273. The data were phase cali- 



10-20 pc scales. Particular attention is devoted to explor- 
ing the influence that counter-rotating instabilities may 
have in driving gas inflow. We also compare in detail the 
CO maps with other gaseous/stellar tracers in order to 
obtain a global picture of the gas response to the stellar 
potential in NGC4826. Finally, we study the distribution 
of star formation in the nucleus, using HST archive broad- 
band and narrow-band images. 

2. Observations 

2.1. CO observations 

Observations of the nuclear region of NGC 4826 were car- 
ried out with the IRAM interferometer between December 
2000 and January 2002, using the BCD set of config- 



S. Garci'a-Burillo et al.: Molecular Gas in NUclei of GAlaxies (NUGA) 



5 




-P I ! ! ' ! ! ' ! ' i ! 1— ' 

D- -* 

Fig. 3. The VLA radio-continuum map at 6 cm from 
Turner and Ho (1994), in grey logarithmic scale from 0.2 to 
12mJy beam -1 . We highlight the position of the dynami- 
cal center derived from CO (star marker) which coincides 
with a strong non-thermal source. 



brated in the antenna-based mode. The flux of the primary 
calibrators was determined from IRAM measurements and 
taken as an input to derive the absolute flux density scale 
in our map; the latter is estimated to be accurate to 10%. 
The bandpass calibration was carried out using 3C273 and 
is accurate to better than 5%. 

The point source sensitivities derived from emission- 
free channels of 10 km s -1 width are 4.0 mJy beam -1 in 
12 CO(1-0) and 5.7mJy beam -1 in 12 CO(2-l). The im- 
age reconstruction was done using standard IRAM/GAG 
software (Guilloteau & Lucas 2000). Unless explicitly 
stated, we used natural weighting and no taper to gen- 
erate the 1-0 line maps with a field of view of 75" and 
0.25" sampling; the corresponding synthesized beam is 
3.0" x 2.3", PA=44°. We used uniform weighting to gen- 
erate 2-1 maps with a field of view of 43" and 0.20" 
sampling; this enables us to achieve a spatial resolu- 
tion <1" (1.1" x 0.8", PA=34°). No 3mm (1mm) con- 
tinuum was detected towards NGC4826, down to an rms 
noise level of 0.38 mJy beam -1 (0.67 mJy beam -1 ) in a 
564MHz-wide band centered on 113.6GHz (231.7GHz). 
The conversion factors between Jybeam -1 and K are 
13 K Jy -1 beam at 115 GHz, and 25 K Jy -1 beam at 
230 GHz. By default, all velocities are referred to v Q and 
(Aa, AS) offsets are relative to the phase tracking cen- 
ter. Except for channel maps, all displayed maps are cor- 
rected for primary beam attenuation. We will assume 
a distance to NGC4826 of D=4.1x(75H -1 ) Mpc (Tully 
1988); H = 75kms -1 Mpc -1 implies l"=20pc. We will 
assume that the inclination angle of NGC4826's disk is 
i=Q0° (Rubin 1994); we derive here a position angle of 



Pj4=112±8°, in rough agreement with previous determi- 
nations (e.g., PA=120±5°; Rix et al. 1995). 

2.2. Optical and near-infrared observations 

We acquired from the HST archive three broad-band im- 
ages of NGC 4826, including WFPC2 (F450W (~ B band) 
and F814W (~ / band)) and NICMOS (NIC3: F160W 
(~ H band)); we also obtained a Paa image (F187N) 
as this galaxy was included in the survey of Boker et al. 
(1999). The optical images were combined with elimina- 
tion of cosmic rays (crreject), and calibrated as described 
in Holtzman et al. (1995). Sky values were assumed to be 
zero. This was checked to be in general a good assump- 
tion with related NICMOS data (Hunt & Malkan 2003 in 
prep), and gives H-band magnitudes that agree to within 
0.1 mag with ground-based data. The NICMOS images 
were re-reduced using the best calibration files and the van 
der Marel algorithm was used to remove the "pedestal" 
effect (see Boker et al. 1999). The continuum was sub- 
tracted from the Paa image according to the precepts ex- 
plained in Boker et al. (1999), but only for the brightest 
pixel values (we excluded roughly half the pixels in the im- 
age). It was necessary to adjust the slope of F160W-F187N 
such that the nucleus remained > everywhere (see also 
Garci'a-Burillo et al. 2000) . The fields-of-view are 2 . 7' x 2 . 7' 
and 51"x51", respectively, for the WFPC2 and NICMOS 
images. 

Optical and near-infrared surface-brightness profiles of 
NGC 4826 were measured along the minor axis (PA = 22°) 
with a 3" width, and used to derive color profiles after 
rebinning to a 0.1 " pixel scale. These will be discussed in 
Sect.S 

3. Results 

3.1. Position of the AG N 

Fig. lll2l show the velocity-channel maps of 12 CO(1-0) and 
12 CO(2-l) emission in the central region of NGC 4826. 
The typical spider web diagram visible in the maps is the 
signature of a spatially resolved rotating disk. The east- 
ern (western) side of the CO disk is blue (red)-shifted with 
respect to the reference velocity v Q . The velocity field in- 
dicates that molecular gas co-rotates with the stars (Rix 
et al. 1995) and the HI gas (Braun et al. 1992) to a ra- 
dius of 700 pc. We find no evidence of counter-rotating 
molecular gas in the inner disk. The best fits for the 
center of rotation and systemic velocity are (Aa,A(S)=(- 
3.3", -0.8") and v^ s s fl =v +8 km s -1 =416 km s -1 , re- 
spectively. Within the errors, the derived dynamical cen- 
ter coincides with the position where CO lines reach 
their maximum widths (^150-200 km s -1 ). The CO dy- 
namical center lies very close (±0.2") to a strong non- 
thermal continuum source present in the 6cm map of 
Turner & Ho (1994) (their source number 5; see also 
Fig. El in this work). Furthermore, the dynamical cen- 
ter coincides, within the errors, with the position of 
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a X-ray source detected in the high-resolution (0.5") 
CHANDRA image of the nucleus of NGC 4826 (source 
S3; PLGordon Garmire, id. 411). Therefore, we ascribe 
the position a J20 oo = 12 h 56 m 43.64 s , 6j 20 oo = 21°40'59.3" 
to both the dynamical center and the LINER nucleus 
in NGC 4826. The systemic velocity derived from CO 
[v sys (LSR)=4l6±5 km s~ 1 =v sys (HEL)=4Q8±5 km s" 1 ] 
agrees satisfactorily with the value determined from HI 
[v sys (HEL)=A08±8 km s" 1 from RC3]. 

3.2. The NGC 4826 s molecular gas disk 

3.2.1. Molecular gas masses 

The velocity-integrated 12 CO(1-0) flux within the 42" pri- 
mary beam field corrected for primary beam attenuation is 
Sco=l-58x 10 3 Jy km s _1 . Assuming a CO-to-H2 conver- 
sion factor X=N(H 2 )/I co =2.2xl0 20 cm- 2 K" 1 km" 1 s 
(Solomon & Barrett 1991), the total H 2 mass derived 
from the interferometer map is M(H2)=2.3x 10 8 Mq. 
Including the mass of helium, the total molec- 
ular gas mass in the 12 CO(1-0) Bure field is 
M 9as =M(H 2 +He) = 1.36xM(H 2 )=3.1xlO 8 M , in good 
agreement with the estimate of 3x10 8 Mq by Sakamoto 
et al. (1999). We estimate that the PdBI map recovers 
close to 100% of the single-dish 12 CO(1-0) flux measured 
by Casoli & Germ (1993) with the IRAM 30m. 

The overall distribution of molecular gas is best seen 
in Fig. 21 which shows the velocity-integrated intensity 
12 CO maps. These were derived by integrating channels 
from v= -200 to 200 km s" 1 . Although the 12 CO(1-0) 
map probably lacks the sensitivity to detect gas emission 
beyond r-35", we find little evidence that NGC4826's 
molecular disk extends significantly outside the 42" pri- 
mary beam (see Fig.Q. This confirms the finding of Casoli 
& Gerin (1993) (see also Sakamoto et al. 1999) who re- 
ported that molecular gas is confined to the inner 1.6 kpc 
of NGC 4826. While our 12 CO(2-l) map lacks sensitivity 
outside the central 500 pc of the disk, it gives a sharp im- 
age of the molecular gas distribution in the vicinity of the 
AGN (Fig. 4b). 

3.2.2. The imprint of m=l perturbations in the 
central kiloparcsec of NGC 4826 

Nearly ~15% of the total molecular gas mass is located 
in a circumnuclear disk (hereafter CND) around the AGN 
(M^ D =3.4x 10 7 M o ). This disk has a ~80 pc average ra- 
dius and its boundaries are here defined by the eighth in- 
tensity contour of Fig. 4a. The small scale structure of the 
CND is fully resolved in the 2-1 map (Fig. 4b). At scales 
of ~tens of pc, the CND has a lopsided morphology: there 
is a pronounced 20-60 pc offset between the AGN locus 
and a ridge of 12 CO(2-l) emission which roughly extends 
from the NE (-3", 0.5") to the W side (-6.5", 1") of the 
CND (see Fig. 4b). A similar offset between the center 
of NGC 4826 and the centroid of dense molecular gas was 
identified by Heifer & Blitz (1997). Their HCN(l-O) map 



revealed a strong emission peak at ~(-1.5",-1.5"), i.e., 
close to the ridge seen in 12 CO(2-l), offering supporting 
evidence that the distribution of molecular gas in the in- 
ner region of the CND (R<60 pc) is strongly lopsided. In 
the outer boundary of the CND (R~60-80 pc) two weak 
winding spiral arcs can be tentatively identified in the 2-1 
map, suggesting that secondary perturbations are at play. 

Outside the CND, molecular gas emission is also de- 
tected along two spiral arms which are located at notice- 
ably different radii in the disk. This indicates that also 
at these scales (100-700 pc) the distribution of molecular 
gas is asymmetrical. An inner spiral arm (hereafter Arm 
I), identified in both the 1-0 and the 2-1 maps (Fig. 
extends south of the CND. The CND and Arm I join at 
~(2.5",-3"). An outer spiral arm (hereafter Arm II), vis- 
ible in the 1-0 map (Fig. 4a), extends north of the CND, 
starting from the western side of the disk [~(-10", 5")]. 
Arm II ends apparently on the SE side of the disk [at 
~(22",-15")]. However, CO emission is detected south- 
ward of the CND at declination offsets ~-20" with a sig- 
nificant lOcr level. Together with Arm II, this S arc would 
delineate an off-centered ring. 

To highlight the geometry of Arms I and II, we have 
constructed the 12 CO(1-0) and 12 CO(2-l) peak bright- 
ness intensity maps (sensitive to contrasted structures) 
and deprojected them into the plane of the galaxy. Results 
are displayed in Fig. [5] which represent the deprojected 
brightness maps in polar grid coordinates [Log e (r),$]. 
Variable r is the galaxy deprojected radius in " and $ 
is the azimuthal angle in degrees (see Fig. for details). 
In this representation, a spiral logarithmic feature would 
appear as a straight line with a non-zero slope. 

As shown in Fig. [S] Arm I develops between 
2.4<Log e (r)<2.8 (ll"<r<16") and 160°<$<280°. Arm 
I shows only a 2 x 7r-periodicity along $. In contrast to 
two-arm (m—2) spiral arms, one-arm (m=l) spirals are 
not 7r-periodic in azimuth. We can derive an average pitch 
angle (p) of the spiral from tan(p)=-ALog e (r)/A$ along 
the arm; p is the angle between the spiral and a circle 
at a given radius, measured counter-clockwise. This gives 
p ^167° for Arm I. The dust lane seen in optical images 
of NGC 4826 is a projected foreground layer situated on 
the northern side of the galaxy, identifying this as the 
nearer side. According to the observed sense of rotation 
of the gas, we therefore conclude that Arm I is trailing 
with respect to the gas and the stellar flows. Interestingly, 
this is in agreement with the morphological analysis by 
Walterbos et al. (1994), who also deduced that the inner 
spiral structure must be trailing inside R^2.7 kpc. 

Arm II is seen in the outer disk in the range 
2.8<Log e (r)<3.5 (16"<r<33") and 65°<$<190°. Arm II 
shows a low order 2 x 7r-periodicity along $, characteristic 
of m=l instabilities. A similar analysis as done for Arm I 
indicates that the pitch angle of Arm II changes midway 
from p ~146° (trailing) to p ~10° (leading). This sudden 
change of orientation might suggest that at this location 
(R~600 pc) two spiral-like perturbations are meeting in 
the disk. 
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In summary, the distribution of molecular gas in the 
inner 1 kpc of NGC 4826 shows the prevalence of different 
types of asymmetrical m—l perturbations in the gas disk. 
Besides two one-arm trailing spirals-Arms I and II-which 
develop in the outer region (radii ~100-700 pc), we have 
detected a lopsided CND. This indicates that m=l pertur- 
bations extend to radial distances as small as ^20-60 pc 
from the AGN. As shown below, the gas kinematics in the 
inner 1 kpc reveal the presence of non-circular motions 
which are related to the various detected m=l perturba- 
tions. A detailed analysis of these is developed in Sect. 
13.3.21 to give an insight into the origin of m=l perturba- 
tions in NGC 4826. 

3.3. Kinematics of NGC 4826's molecular gas disk 

3.3.1. The rotation curve and the dynamical mass 

We show in Fig. El the mean velocity field maps obtained 
from the 12 CO(1-0) and 12 CO(2-l) data in the nucleus of 
NGC 4826 (derived with 2-cr clipping). The kinematics of 
molecular gas at radii R<700 pc are consistent with those 
of a disk in direct rotation with respect to the stars. 

A CO rotation curve (v ro t) has been derived from 
position- velocity (p-v) diagrams taken along the kinemat- 
ical major axis of NGC 4826 (Fig.Q). We identify the ma- 
jor axis position angle as PA=112±8°; this value yields a 
larger line-of-sight velocity gradient within R=700 pc than 
does a p-v cut through the dynamical center at any other 
angle. Our determination roughly agrees with previous 
findings based on Ha and stellar kinematics (PA=120±5°; 
Rix et al 1995). We have calculated the terminal veloci- 
ties, and from these, the rotation curve, by fitting gaussian 
profiles to the spectra across the major axis. The veloc- 
ity centroids, corrected for inclination i=60°, give v rot for 
each offset along the major axis. Rotation curves derived 
from both sides of the major axis do not differ signifi- 
cantly within the errors. Therefore, we have derived v rot 
as a function of radius by averaging values from the W and 
E sides of the major axis. Data for both lines were noise- 
weighted averaged using a radial beaming of Ar=0.5". 

The resultant v rot (Fig. |HJ) shows a steep increase to 
the edge of the CND (R~70-80 pc), where v rot ~120- 
125 km s _1 . This trend is followed by a gradual increase 
up to a radius of R~350 pc, after which v rot remains 
flat at ^190 km s -1 . The rotation velocity of molecular 
gas follows the stellar and the ionized gas motions inside 
R-700 pc (see Fig. 3 and 6 of Rix et al 1995). 

The total mass inside the CND can be inferred using 
M(R)=CxRxv 2 ot /G, where G is the constant of gravita- 
tion, M(R) is the mass inside a sphere of radius R, and 
C is constant varying between 0.6 and 1, depending on 
the disk mass model assumed (Lequeux 1983). If we take 
a value of C=0.8 intermediate between the values appro- 
priate for spherical (1) and flat disk (0.6) distributions, 
M(R=70 pc)~2.3xlO 8 M . This implies that the molecu- 
lar gas mass fraction inside the CND is ~15%, decreasing 
to ~5% within R=700 pc. 



With the value estimated for the rotation curve at a 
radius R~10pc (v rot ~58 km s _1 ) we can set an upper 
limit of ~8x10 6 Mq for the mass of the putative super- 
massive black hole in the nucleus of NGC 4826. 

3.3.2. Streaming motions and m=l modes 

Although mainly characterized by circular rotation, the 
gas kinematics are perturbed by streaming motions. 
Isovelocities in Fig. display a wavy pattern at the pas- 
sage of both Arms I and II. Deviations from circular mo- 
tion can be identified in Arm I as a systematic kink in 
the isovelocities which appear redshifted; as expected, the 
shift is most pronounced near the crossing of the galaxy 
minor axis where only radial motions have non-zero pro- 
jection (see Fig. EJl . We find a gradient across the minor 
axis from redshifted velocities in Arm I to blueshifted ve- 
locities in the adjacent interarm region located at larger 
radii (see Fig. 6a and Fig. [SJ. Deprojected into the galaxy 
plane, the radial velocities measured on Arm I indicate 
local outflow motions. In contrast, non-circular motions 
along Arm II indicate that gas velocities are systematically 
redshifted, a local signature of inflow motions (Fig. [HI- 

Departures from pure rotation, related to the lopsided 
instability, also characterize the CND kinematics. The av- 
erage gas velocities derived from 2-1 data seem slightly 
redshifted (blueshifted) on the southern (northern) cross- 
ing of the lopsided feature along the minor axis (Fig. |5J). 
Similarly to Arm I, these motions can be interpreted as a 
local signature of outflow. However radial streaming mo- 
tions are neatly reversed in the outer boundary of the 
CND (R~60-80 pc) where two secondary winding spiral 
arcs are identified in the 2-1 map (see Fig.[!|J). 

The signatures of the streaming motions expected for 
trailing/leading waves can be analyzed in the framework of 
the linear density wave theory (Shu et al. 1973). The sign 
of velocity perturbations changes close to the minor axis 
when crossing from inside to outside the corotation reso- 
nance. A comparison with models may give some insight 
into the nature of the density waves that may account for 
the streaming motions observed in CO. Garci'a-Burillo et 
al. (2000) applied the general formalism to the case of one- 
arm (m=l) spirals for two extreme values of the pattern 
speed (fip) illustrating the slow and the fast mode solu- 
tions. In the slow (fast) solution we are mostly inside (out- 
side) corotation of the modes. In the case of NGC 4826, 
we have explored a fast solution by adopting a pattern 
speed large enough to assure corotation lies well inside 
the disk of the 'model' galaxy (f2 p =800 kms -1 kpc -1 ; 
this implies Rco/J=200pc for a rotation curve similar to 
v rot ). In contrast, a slow mode is here characterized by 
fl p =0 kms^kpc" 1 , i.e., a stationary wave where corota- 
tion is pushed outside the disk. 

We represent in Fig. ^] the radial velocity perturba- 
tions, projected into the plane of the sky (v pert ), for the 
fast and slow one-arm trailing solutions. We purposely 
chose a disk geometry and spiral parameters qualitatively 
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similar to those of NGC4826, apart from a lower inclina- 
tion angle adopted here to give a more detailed picture of 
the velocity field in the model. Moreover, we can reason- 
ably assume that the gas response (peak gas density) is 
close to the spiral potential minima in all cases. A compar- 
ison of the two model solutions with the CO observations 
indicates that the streaming motions measured on Arm 
I and, also, on the inner lopsided instability of the CND 
match qualitatively the solution of the fast trailing m=l 
mode (Fig. I10J1 : when the spiral arm is outside corota- 
tion, velocities are redshifted (blueshifted) on the south- 
ern (northern) side of the galaxy minor axis. In contrast, 
Arm II better matches the solution of the slow trailing 
mode (Fig. EH. 

A paradoxical consequence is that the inner m—l per- 
turbations, represented by Arm I and the lopsided CND 
instability, would not favor gas infall: behaving as fast 
modes, the gas would gain angular momentum from the 
waves and would migrate outwards. In particular the AGN 
fueling might be temporarily blocked upon the onset of the 
inner lopsided instability. The perturbations identified in 
the outer boundary of the CND would favor gas inflow, 
however. This suggests that inflow and outflow may be 
globally counterbalanced in the CND. It remains to be 
determined, however, if angular momentum removal/gain 
from these instabilities, either stellar or gaseous, can be 
efficient enough to play some role in AGN feeding (see 
Sect.EJ. 

3.4. CO line ratios 

To probe the physical conditions of the molecular gas in 
NGC4826, we have studied the variation of the 12 CO(2- 
1)/ 12 CO(1-0) ratio in areas with significant emission lev- 
els in both lines (brightness temperatures >5er). The 2-1 
map was first degraded to the spatial resolution of the 1- 
map within the 21" primary beam to assure that the 
two lines sample identical regions; both maps were also 
corrected for primary beam attenuation. 

The derived 12 CO(2-1)/ 12 CO(1-0) ratio is roughly 
constant over most of the disk, with a mean value 
0.65±0.05. This ratio is close to the canonical value ob- 
served in other spiral disks and is typical of moderately 
dense optically thick molecular clouds (see Garcfa-Burillo 
et al. 1993). The only exception to this is in the vicinity of 
the AGN (R<50-100pc), where the measured ratio is close 
to 1. The highest values (~1.1) are observed toward the 
lopsided structure in the CND which coincides with the 
peak of emission detected in HCN (Heifer & Blitz 1997). 
Ratios are close to 0.8-0.9 in the strongest clumps along 
the spiral ridge of Arm I. 

4. Star formation in NGC4826 

4.1. Distribution of star formation 

The distribution of star formation in NGC4826's disk is 
strongly asymmetrical, as seen in the HST Paa image 



(see Fig. Illfl . The asymmetries follow with similar spa- 
tial scales the various m=l instabilities identified in the 
CO maps. Interestingly, the overall star formation pat- 
tern shows a strong N/S asymmetry. While the 12 CO(1-0) 
emission on Arm II is associated with a strong ionized gas 
emission, to the southern far side in Arm I, where extinc- 
tion is smaller, Paa is significantly weaker. This is contrary 
to what would be expected if the reported N/S asymme- 
try was due mainly to heavy extinction by the conspicuous 
northern dust lane. Furthermore, the star formation pat- 
tern is also very asymmetrical in the CND fFig.lll[). We il- 
lustrate this in Fig. El which compares the Paa morphol- 
ogy with the mostly thermal radio continuum emission at 
2cm in the nucleus of NGC 4826 (from Turner & Ho 1994). 
A strong maximum is visible in both maps, which indicates 
that the star formation peaks neither on the AGN nor on 
the lopsided instability seen in the 12 CO(2-l) map, but 
rather 20 pc northward. 



4.2. Extinction in the "Evil Eye" 

The conspicuous dust lane identified in the optical pic- 
tures of NGC 4826 is expected to significantly screen any 
emission to the N. With the assumption of intrinsic ax- 
isymmetry for the stellar disk/bulge, Witt et al. (1994) 
determined the optical depth in the "Evil Eye" . Following 
Witt et al. (1994), we have re-examined the issue of ex- 
tinction by extracting surface brightness profiles along the 
minor axis in F450W, F814W, and F160W. The relation- 
ship of this "cut" to the galaxy morphology is shown in 
Fig. El the profiles, averaged over 3" perpendicular to the 
direction of the cut, are shown in Fig. El The figure shows 
clearly the effects of extinction towards positive (northern) 
radii, and confirms the assumption that the stellar popula- 
tions in the disk of NGC 4826 are rather symmetric, as the 
discrepancy between the northern and southern halves of 
the surface brightness profiles decreases substantially with 
increasing wavelength. However, the ionized gas emission 
is clearly intrinsically asymmetric since it is more promi- 
nent on the near side where there is more extinction. 

We have derived the extinction map using the mean 
colors of the disk, as determined from the outer region of 
the cuts, toward the S, where (broadband) extinction is 
low or zero (see Witt et al. 1994). The B — I mean color of 
the outer (S) region is 1.62, similar to, although 0.16 mag 
bluer than, the equivalent color from de Jong (1996) for 
integrated colors of similar morphological types. With the 
extinction coefficients given in Holtzman et al. (1995) for 
the WFPC2 filters, we then derived the extinction from 
the B — I map. The mean extinction measured in the dust 
lane turns out to be not greater than Ay =1.5 mag, in 
agreement with the findings of Walterbos et al. (1994) and 
Witt et al. (1994). This confirms that the Paa image of 
NGC 4826 is virtually extinction-free, and can be consid- 
ered as a fair unbiased picture of how recent star formation 
proceeds in the central kpc of NGC 4826. 
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Fig. 13. a)(left)The F450W image of the galaxy with the "cut" aperture superposed. The brightness peak is shown 
with a cross in the black (brightest) central region, while a bright HII region is marked with a cross at a radial distance 
of 17". b) (right) The surface brightness profile along the minor axis, together with the F450W-F814W ~ B — I and 
F450W-F160W ~ B — H colors. The center and HII region at r ~ 17" shown by crosses in a) (left) are denoted 
here by vertical dotted lines. The F160W profile is shown by a dotted line, F814W by a solid line, F450W by a dashed 
line, and CO(l-O) by a heavy line (from top to bottom). Y axis units are magnitudes for all profiles, 'including' CO 
(in logarithmic scale with an arbitrary offset). In the top panel (for all curves except that for F450W), we plot the 
reflection of the southern half of the surface brightness profile underneath the northern half of the profile, showing 
the asymmetry induced by the dust lane. The horizontal dashed lines in the lower panels show the intrinsic color as 
estimated from the southern region, i.e., at negative offsets along the cuts. 



The average extinction in the nucleus of NGC 4826 can 
also be estimated from the 12 CO(1-0) integrated intensity 
map of Fig. Assuming a CO-to-Eb conversion factor 
A=N(H 2 )/Ic7O=2.2xl0 20 cm" 2 K" 1 km" 1 s (Solomon & 
Barrett 1991) and the widely used gas-to-dust ratio of 
N(H 2 )/Ay = 1.0 10 21 cm" 2 mag" 1 (Bohlin et al. 1978), 
we obtain a mean Ay ~ 10 mag within the 42" 1-0 pri- 
mary beam. This value is a factor of 6 higher than that 
derived from B — I color. However, for several reasons 
the CO-based Ay value likely overestimates the dust col- 
umn density that can absorb efficiently the stellar light. 
First, when this comparison is made, we implicitly assume 
a screen geometry for the H 2 slab with unit filling factor. 
Second, there is evidence that the CO-to-ifa column den- 
sity ratio X in central regions of galaxies may be lower 
than the standard Galactic value (see Garcfa-Burillo et 
al. 2000 and references therein). The extinction 'discrep- 
ancy' is even more severe at small scales. Fig.ll4lcompares 
the extinction map, estimated from B — I, with the CO 
maps. Confirming the impression given in Fig. 1131 it is 
evident that molecular gas column density is not tightly 
correlated with the B — /-based Ay (see also Block et al. 
1994 and Sakamoto et al. 1999). Moreover, the clumpy 
structure of the interstellar medium makes difficult any 



comparison between Ay values derived using data taken 
at very different spatial resolutions. 

4.3. The star formation rate 

We have measured the Paa flux within the entire NICMOS 
field-of-view for NGC 4826, summing > 3er above sky, and 
found F(Paa) = 7.4 x lO -13 erg/s/cm 2 . This gives a lumi- 
nosity L(Paa) = 1.5 x 10 39 erg/s, or a massive star forma- 
tion rate (SFR) of 0.14M Q /yr, after correcting for a mean 
extinction Ay of 1.5 mag (this work). This is slightly lower 
than, although comparable to, the SFR of 0.22 M©/yr es- 
timated by Braun et al. (1994) from an Ha image with a 
larger field-of-view. While this indicates rather vigorous 
star formation activity, it is much lower than in true nu- 
clear starburst galaxies such as M82, which have SFRs 
more than 10-50 times higher (O'Connell et al. 1995). 

The number of ionizing photons can also be measured 
from the ionized gas emission, and from this we can de- 
rive the number of equivalent O stars necessary to ion- 
ize the gas. Using the coefficients given in Ostcrbrock 
(1989), we find that the total Pact luminosity (not cor- 
rected for extinction) corresponds to 9.3 x 10 51 ionizing 
photons. Taking an 07V star as representative, and as- 
suming the number of Lyman continuum photons from 
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such a star to be 1 x 10 49 (Leitherer 1990), this corre- 
sponds to roughly 930 07V stars. Our estimates of ion- 
izing photons are consistent with those of Pierini et al. 
(2002), if we consider the total area of the Pact image 
compared to that of their spectroscopic slit. Because ex- 
tinction may be an issue here, it is important to compare 
the number of ionizing photons derived from Paa to that 
derived by Turner & Ho (1994) from the radio contin- 
uum. They find a factor of 3 fewer ionizing photons, but 
20 times more massive stars, differences which almost cer- 
tainly are due to the different normalizations. As noted 
by Turner & Ho (1994), the HII regions in NGC4826 are 
one order of magnitude more energetic than the largest 
star-forming complexes in the Milky Way. In any case, 
the rough agreement between the estimates obtained from 
Paa and from radio continuum maps is an independent 
confirmation that Paa suffers from virtually no extinction 
(see also discussion in Sect. 14.21 . 

4.4. Star formation and compact star clusters in 
NGC4826 

The comparison of the Paa morphology with that of the 
F450W and the B — I color map reveals an interesting fea- 
ture. The very blue compact source that emerges in the 
color cut at 17" lies within a ridge of HII regions or com- 
pact star clusters, visible in Fig. El just inside the edge of 
the "Evil Eye" . Surprisingly, this ridge of star formation 
does not coincide with the ionized gas emission ridge seen 
in Paa (Fig. Illfl . The Paa ridge is positioned roughly 4" 
closer (80 pc) to the nucleus than the arc of clusters, while 
the (very blue) 17" feature seems to contain little if any 
ionized gas. This result is consistent with the radial vari- 
ations of ionized gas emission along the slit of Pierini et 
al. (2002), who found the peak of the ionizing photon flux 
at ~ 14 — 15"; at 17" the ionizing flux is approximately 
30% of its peak value. These authors also found the hard- 
ness of the ionizing ratio to vary radially, with an abrupt 
softening at 17", roughly the same location where we find 
the blue compact source in the color cuts. 

To better understand the pattern of star formation in 
NGC4826, we have measured F450W and F814W mag- 
nitudes and (FWHM) diameters of compact sources in 
the ridge just inside the edge of the dust lane (see extra- 
nuclear cross in Fig. I13fl . All of the 25 sources measured 
appear to be star clusters, since their absolute F814W ~ I 
magnitudes range from —7.4 to —8.9 (without correction 
for extinction) and their FWHMs from 4 to 29 pc. If we 
convert to V assuming B — V ~ 0.7 and correct for 1.5 
mag of visual extinction, we find that the absolute V mag- 
nitudes (—8.5 to —9.9) are consistent with those of Super 
Star Clusters (SSCs) if they are old enough (in which case 
the correction to the fiducial age of lOMyr would be > 
1 mag) . 

If we infer the ages of these clusters from their B — I 
color, which ranges from 0.65 to 2, we would derive a min- 
imum age of 40 Myr (Leitherer et al. 1999) in the absence 



of reddening. This would imply a minimum luminosity 
correction of ~ 1.5 mag (Whitmore et al. 1997), which 
would bring the luminosity of most of these compact star 
clusters into the SSC range of My < —10.5 at a fiducial 
age of 10 Myr (Billet et al. 2002). These clusters are not 
extremely compact, although this also could be related 
to an older age. In any case, the clusters must be older 
than 10 Myr, given the conspicuous absence of ionized 
gas emission in their immediate vicinity. The formation 
scenario for compact star clusters or SSCs, although not 
completely understood, appears to involve mergers (Bekki 
& Couch 2001) and/or high-pressure environments (e.g., 
Billet et al. 2002). 

The ionized and the molecular gas profiles along the 
minor axis of NGC 4826 share a common feature: the disks 
are abruptly truncated beyond a ridge of old (>10Myr) 
compact star clusters. This pattern suggests a parallel 
pattern in the age of the star formation in the disk. 
Apparently stars are still vigorously forming elsewehre in- 
side R^700pc where there is a significant molecular gas 
reservoir. 

5. Counter-rotation and m=l modes 

NGC 4826 is a prototype among the large variety of 
counter-rotating disk galaxies. Counter-rotation can be 
purely stellar, involving stars versus gas, or character- 
ized by the co-existence of two decoupled gas disks as 
seen in NGC 4826 (see compilations by Bertola & Corsini 
1999, and Bettoni et al. 2001). The existence of counter- 
rotating components in disks suggests the active role that 
interactions, minor mergers and mass accretion processes 
may have in driving the evolution of galaxies. In hierar- 
chical merging scenarios which explain how galaxy disks 
are built up, counter-rotating galaxies may naturally re- 
sult from merger and/or accretion events. Partly due to 
the small number of observations, the prevalence of decou- 
pled components in disk galaxies may have been underes- 
timated thus far. Nevertheless, it is important to investi- 
gate which are the main gravitational instabilities induced 
by counter-rotation, and to analyze their influence on the 
removal of angular momentum from the gas. Analytical 
studies by Lovelace et al. (1997) indicate that the main dy- 
namical instabilities linked with two-stream flows in galac- 
tic disks are one-arm spirals, i.e., m—1 modes. 

Garcfa-Burillo et al. (1998, 2000) mapped at high- 
resolution the massive counter-rotating molecular disks of 
the spirals NGC 3593 and NGC 3626, searching for the sig- 
natures of counter-rotating instabilities. The maps of both 
galaxies show compact molecular gas disks with mixtures 
of m=l and m=2 perturbations. With the help of self- 
consistent numerical simulations adapted for NGC 3593, 
Garcfa-Burillo et al (2000) studied the development of 
disk instabilities. In their study, they found that counter- 
rotation drives rapidly evolving perturbations. The disk 
first develops stationary (O p ~0 kms -1 kpc -1 ) waves, 
leading with respect to the gas flow, followed later by 
a mixture of m=l and m=2 slow (trailing) waves. The 
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nature of the instabilities is seen to depend critically on 
the assumed halo/disk mass ratio (see also Comins et al. 
1997). Similarly to the aforementioned cases, NGC4826 
has a very compact molecular disk with signatures of m=l 
perturbations at various scales. Rix et al. (1995) analyzed 
the stellar velocity distributions and estimated that a siz- 
able fraction (10%-30%) of the stars in NGC 4826's inner 
disk (radii <50") could be in counter-rotation. This decou- 
pled component could drive the two-stream flow instabil- 
ities we actually observe in the molecular disk. However, 
the reported fraction of counter-rotating stars is still com- 
patible with the high-velocity dispersion of this kinemati- 
cally hot region, and thus cannot be taken as firm evidence 
for counter-rotation in the stellar component. 

In summary, the only true counter-rotating component 
in NGC 4826 may be the outer HI disk. The instabili- 
ties observed in the inner molecular disk are therefore not 
linked with counter-rotation: m=l perturbations are de- 
tected in a region where two-stream flow is not detected. 
Additional evidence for this comes from the analysis of 
streaming motions made in the frame of density wave 
models (see Sect. l3.3!2*)l . The CND lopsided instability and 
Arm I behave as fast modes, i.e., if the inner m=l per- 
turbations are wave modes they seem to have developed 
outside corotation. This contrasts with the low pattern 
speeds of m=l modes which are known to be generated 
by the two-stream flow instability (Lovelace et al. 1997; 
Comins et al. 1997; Garci'a-Burillo et al. 2000). 

6. Gravitational instabilities and AGN fueling 

6.1. Gravitational torques in NGC4826 

The objective of this section is to study if gravitational 
torques, derived from a fair representation of the stellar 
potential in the inner disk of NGC 4826, can account for 
the gas kinematics derived from CO. Furthermore, we ex- 
amine the efficiency of gravitational torques exerted on 
the gas in the context of AGN fueling. 

NIR maps can be used to derive the distribution of 
old stars as they are less affected by dust extinction or 
by stellar population biases. Here we have used a high- 
resolution near-infrared image of NGC 4826 taken with 
HST-NICMOS to derive the stellar potential within the in- 
ner 1 kpc of the disk. The HST field-of-view on NGC 4826 
is 1024 pc and is sampled with a grid of 256x256 pixels 
of 0.203" spatial resolution. The image of the galaxy was 
completed along the 3rd vertical dimension by assuming 
an isothermal plane with a scale height constant with ra- 
dius; this scale height is equal to l/12th of the radial scale- 
length, i.e.,~160 pc. The gravitational potential is derived 
by a Fourier transform method on the grid from the NIR 
image. We also assumed a constant mass-to-light (M/L) 
ratio; the value of M/L is obtained by fitting the observed 
rotation curve. 

The potential is decomposed as 



we define the strength of the m-Fouricr component, 
Qm{R) = m& m /R\Fo(R)\ and its global strength over the 
disk as max#Q m (i?) (e. g. Combes & Sanders 1981). 

The strength of the total potential or maximal torque 
over the whole disk is defined by 
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where F?p ax (R) represents the maximum amplitude of the 
tangential force and Fq(R) is the mean axisymmetric ra- 
dial force, inferred from the m=0 component of the grav- 
itational potential. 

The strengths of the main m components for R<520 pc 
are plotted in Fig. together with the radial-variation 
of the corresponding phases, </> m . It is interesting to notice 
that m—l components are stronger than m=2 components 
in the stellar potential inside this region, especially for 
R<50 pc. This indicates that there are also asymmetrical 
perturbations in the stellar disk. 

Once the potential is obtained all over the grid, we 
deduce the forces {F x and F y ) by derivation in each pixel, 
and can compute the average torque exerted on the gas. 
The surface density E of the gas is assumed proportional 
to the CO emission, either 12 CO(1-0) or 12 CO(2-l) (both 
maps are used). We resample the CO maps at the same 
pixel size of 0.203", and deproject them with the same 
geometrical angles as before. On each pixel, the torque T 
is computed by: 

T = Z(x.F y - y.F x ) 

We then compute the average over each radius of the 
torque per unit mass, i.e.: 
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The mean torques calculated using the 12 CO(1-0) map 
are plotted in Fig. ^1 in units of (km/s) 2 . The results 
based on the 12 CO(1-0) and 12 CO(2-l) maps are consis- 
tent. Most remarkably, oscillations between negative and 
positive values for t(R) are roughly consistent with the di- 
agnostic based on the change of sign of streaming motions 
for Arm I (related to a fast trailing wave, hence causing 
angular momentum gain) and Arm II (related to a slow 
trailing wave, hence causing angular momentum loss) (see 
Sect. l3.3^|l . The corresponding stellar torques in the CND 
are marginally positive very close to the AGN (R<50 pc), 
suggesting that stellar perturbations contribute little to 
AGN feeding. Again, this fits 'qualitatively' the stream- 
ing motions measured in CO, which would rather sug- 
gest outflow to some degree. The stellar torques change to 
negative values in the outer boundary of the CND; this 
roughly agrees also with the CO-based diagnostic. 

While the radial variation of the stellar torques seem to 
account qualitatively for the changing signature of stream- 
ing motions, the maximum value of the mean torque is ex- 
ceedingly small, however: ^50(fcm/s) 2 . This implies that 
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the typical time-scale for the gas to lose its angular mo- 
mentum for R~200pc is of about 1 Gyr. This is much 
longer than the gas dynamical time-scale at this radius (of 
the order of 2 Myr). Gravitational torques exerted on the 
gas due to the stellar potential are 'globally' quite weak. 
Weighted by the density, the total torque inside the evalu- 
ated region of the disk (R~520 pc) is only marginally neg- 
ative: the combination of all stellar perturbations seems 
to make the gas lose its angular momentum rather ineffi- 
ciently. 

We can ask whether our conclusions might be partly 
biased by the existence of residual extinction in the H- 
band image originally used to derive the gravitational 
potential. However, several arguments can be advanced 
suggesting that this bias is negligible. Although we es- 
timate that extinction in the H band can locally reach 
Ah ~ 0.20 — 0.25 magnitudes in the northern dust lane 
arc (at R ~ 250 — 300 pc), this value is comfortingly low 
elsewhere. In particular, we do not expect that the stellar 
gravitational potential is severely biased inside the CND, 
i.e., at the scales which are critical for evaluating the in- 
fluence of gravitational torques on the AGN feeding. The 
derived values for t(R) inside the CND (at R < 80 pc 
where Ah < 0.02 — 0.03), as well as in Arm I (develop- 
ing on the southern side, where Ah < 0.01), are mostly 
weighted by the contribution of regions that are virtually 
extinction- free (see Fig. 1140 . 

To estimate quantitatively the global influence of ex- 
tinction on the derived torques, we have corrected the NIR 
map by a factor exp(r^f), where tjj = O.lAAy and Ay 
corresponds to the extinction map of Fig. 1141 Confirming 
our expectations, the Fourier analysis of the gravitational 
potential thus derived is only slightly modified, and the 
average torques on the gas are practically unchanged. 
Although a detailed evaluation of residual extinction on 
the 7J-band image would require the use of several high- 
resolution optical and infrared images in order to resolve 
the age-extinction degeneracy (to be addressed in a forth- 
coming publication), we can conclude that the value de- 
rived here for t(R) in the center of NGC 4826 is a robust 
estimate. 

6.2. Gas self-gravity 

Stellar perturbations are weak in the inner disk of 
NGC 4826 (except in the very center of the galaxy: 
R<50 pc). Furthermore, we have found indications that 
they are probably inefficient to drive AGN fueling at 
R<520 pc. Alternatively, self-gravitating gas instabilities 
may have a potential role in AGN feeding. When self- 
gravity of stars and/or gas is strong, galaxy disks are prone 
to gravitational instabilities. The local stability crite- 
rion, first established by Toomre (1964), ensures stability 
against axisymmetric perturbations in a stellar/gaseous 
galaxy disk when the column density of stars/gas (iV* ifl ) 
is smaller than a critical value (N% r g t =Kcr v /%G) , i.e., when 
the so-called Toomre parameter Q*,g=N%™ t /N* ig > 1. 



Here, k is the epicyclic frequency, a v is the velocity disper- 
sion and G is the gravitational constant. Stability against 
non- axisymmetric perturbations, such as spirals and bars, 
is guaranteed when Q is significantly larger than 2. 

The stability criterion for a realistic galaxy disk needs 
to take into account the strong coupling between the stel- 
lar and gaseous components through self-regulation and 
feedback (Combes 2001). Nevertheless, as a first step, it is 
useful to evaluate Q separately for the gas disk. In the 
particular case of AGN host disks, previous extinction 
maps derived from HST color images seem to indicate 
that non self-gravitating gas instabilities (Q ~10-100) may 
be ubiquitous in the inner 100-200 pc of Seyfert galaxies 
(Martini & Pogge 1999) and potential drivers of AGN fu- 
eling through energy dissipation in shocks and turbulent 
motions. This scenario can be fully tested by exploiting 
the information provided by 12 CO maps. Under certain 
assumptions, the 12 CO intensity gives a fair estimate of 
N g . Moreover, the 12 CO kinematics can be used to derive 
f2, k, and a v in the disk. Hence, the derivation of Q maps 
follows straightforwardly. 

Fig. 1 161 shows the variation of Q in the molecular disk 
of NGC 4826 as derived from 12 CO(1-0). Fig. [TB] displays 
also a zoomed-in view of the Q map for the inner disk 
obtained from 12 CO (2-1). Gas column densities N g are 
derived from 1-0 data as explained in Sect. 13. 2. H We have 
used the 2-1/1-0 line ratio maps discussed in Sect. 13.41 to 
calculate N g from 12 CO(2-l). fl and k are derived from 
v rot using standard definitions. Finally, a v has been esti- 
mated from CO second moment maps. The contribution 
to a v of the rotation curve gradient within the lobe(s) is 
found to be significant for r < 2", but negligible elsewhere. 
The velocity dispersion is fairly constant in NGC4826's 
disk and close to 10 km s _1 except for r < 2", where a v 
doubles its value on average. In the following analysis we 
take for simplicity a v = 10 km s . 

With these assumptions, Fig. indicates that Arms 
I and II are self-gravitating perturbations: the two one- 
arm spirals are characterized by low Q values, ranging 
from 1 to 2. There are also indications for the CND lop- 
sided instability to be self-gravitating (Q <2-3), unless 
an unrealistic boost in velocity dispersion (not observed) 
or a much lower CO-to-H2 conversion factor is invoked 
for r < 2". NGC 4826 seems to be a good counterexam- 
ple to the scenario suggested by Martini & Pogge 1999: 
the m=l instabilities identified in NGC 4826 have low Q 
values. This suggests that self-gravity of the gas cannot 
be neglected and probably plays an essential role in the 
maintenance of the m=l-type perturbations described in 
NGC 4826. 

7. Summary and conclusions 

We summarize the main results obtained here as follows: 

— High-resolution 12 CO observations of NGC 4826 show 
that the bulk of its molecular gas lies in a highly struc- 
tured disk of M gas ~3.1x10 8 Mq which ends abruptly 
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at R/^700 pc. The gas disk shows a preponderance of 
asymmetric perturbations (777=1 instabilities). There 
is a lopsided nuclear disk of 40 pc radius and two one- 
arm trailing spirals, which develop at different radii in 
the disk (Arm I: 200-350 pc, Arm II: 350-700 pc). 

— The star formation pattern in the disk of NGC 4826 
is strongly asymmetrical. The scales of the observed 
asymmetries resemble those of the various m=l per- 
turbations revealed in the disk of molecular gas. 
Furthermore, the star- forming ionized/molecular gas 
disk is truncated beyond a ridge of compact star 
clusters identified at R~700 pc. This suggests an 
evolutionary trend for star formation in the disk of 
NGC 4826. Apparently massive star formation has 
ceased beyond R^700 pc, but is still vigorously occur- 
ring inside this radius, fed by the significant molecular 
gas reservoir. 

— Gas kinematics reveal streaming motions related to the 
r7i=l perturbations. A first analysis of these perturba- 
tions suggests that the inner m—l instabilities may 
be fast trailing waves which have developed between 
corotation and an Outer Lindblad Resonance (OLR). 
This would imply that the AGN is probably not be- 
ing generously fueled in the current epoch. An esti- 
mate of the radial variation of the mean gravitational 
torques due to the stellar potential confirms indepen- 
dently that stellar perturbations are inefficient to drive 
AGN fueling. 

— Arms I and II are both self-gravitating perturbations 
characterized by low values of the Toomre Q parame- 
ter (Q=l-2). There are also indications that the CND 
lopsided instability is self-gravitating (Q <2-3). 

The two-stream flow instabilities, expected in counter- 
rotators, cannot explain the ubiquity of m=l pertur- 
bations in the inner disk of the Evil Eye. Mechanisms 
other than counter-rotation have been suggested to trig- 
ger m=l instabilities, such as interactions with compan- 
ions (Weinberg 1994; Lovelace et al. 1999), central poten- 
tials dominated by massive black holes (Miller & Smith 
1992; Taga & lye 1998; Bacon et al. 2001) and response 
to an asymmetric halo (Jog 1997). Although pure 777=1 
modes are uncommon in galaxy disks, kinematic lopsided- 
ness has been reported in 50% of spirals (Richter & Sancisi 
1994). The self-consistent numerical simulations made by 
Junqueira & Combes (1996) showed that above a thresh- 
old central gas concentration, a typical galaxy disk (in- 
cluding stars and gas) is prone to develop one-arm (trail- 
ing) spiral perturbations. The strongest m=l modes in 
their models appear between their corotation and their 
OLRs; in other words, these perturbations are fast. Fast 
m=l trailing modes are decoupled from outer disk pertur- 
bations which only have 1/10 the speed of the inner waves. 
According to Junqueira & Combes (1996), the fast m=l 
modes might reflect the action of the modal amplification 
mechanism described by Shu et al. (1990). The modes are 
excited by a slight off-centering of the stars and gas in the 
galaxy nucleus and develop mainly in the gas disk between 



corotation and OLR. In summary, the inner m—l pertur- 
bations identified in the gas disk of NGC 4826 could be 
tentatively explained by this scenario if the pattern speed 
of the perturbation is tuned to high enough values. As 
shown in Fig. El a value of £l p >1500 km s _1 kpc -1 
would push corotation well inside the disk, accounting for 
the spatial extent of Arm I in NGC 4826. 

7.1. A possible scenario for the evolution of the inner 
kiloparsec of NGC 4826 

While it appears that counter-rotation is not driving 777 = 
1 modes in the inner 1.5 kpc disk of NGC 4826 at present, 
a past gas accretion episode must be invoked to explain 
the decoupling of HI in the outer disk. Furthermore, the 
formation of compact/truncated star forming disks in 
counter-rotators, such as the one we see in NGC 4826, 
could represent the final stage of a process involving large- 
scale collisions between the accreted gas and the primary 
gas of the accretor (Rubin 1994, Thakar et al 1997, Garcia- 
Burillo et al 2000). In general, the two components may 
have opposite angular momenta, implying substantial dis- 
sipation when they mix and, eventually, transformation 
from atomic to molecular phase in shocks. Large amounts 
of gas may fall towards the nucleus and form a circum- 
nuclear gas disk, whose final size will depend on the ra- 
tio of initial angular momenta of the two components. In 
close agreement with this picture, the kinematics of ion- 
ized gas reveal an orderly infall of gas from R=2 kpc to 
800 pc in NGC 4826 (Rubin 1994). Moreover, the mea- 
sured [Nil] /Ha ratio is close to 1 in this transition region, 
strongly suggesting shock excitation (Rubin 1994). In the 
course of this process, massive star formation is triggered 
along the "ridge" and in the nuclear disk. The time-scale 
for gas infall may be short and similar to the dynamical 
time, i.e., ~10 7 ~ 8 yr. This estimate is consistent with the 
age of the stellar clusters inferred from their B — I colors. 

The stellar velocities measured in the nucleus of 
NGC 4826 also imply a recent large mass infall episode. 
Kormcndy (1993) has discussed NGC 4826 as an ex- 
ample of an early-type spiral with an anomalously 
low stellar velocity dispersion for its bulge luminosity 
(<7t,=90±5 km s _1 , later confirmed by Rix et al 1995). The 
expected value of a v for NGC 4826 predicted by the Faber- 
Jackson relation would be ^160 km s _1 . As argued by 
Kormendy (1993), this indicates that the central bright- 
ness is dominated by a cold disky component. The forma- 
tion of this cold central disk, in contrast to standard hot 
bulges, might be related to secular evolutionary processes 
involving large mass accretion. 

Once the gas has settled in the nuclear disk, normal 
secular evolution can proceed. The onset of m=l instabil- 
ities of the type described in NGC 4826 may be a conse- 
quence of secular evolution in disks with large gas masses. 
Detailed numerical simulations to be presented in a forth- 
coming paper (Garci'a-Burillo et al. 2003 in prep) will 
study the onset and evolution of asymmetric modes for 



14 



S. Garci'a-Burillo et al.: Molecular Gas in NUclei of GAlaxies (NUGA) 



a case similar to NGC4826, where the role of self-gravity 
of the gas may be essential. 
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Fig. 4. a) 12 CO(l-0) integrated intensity contours in the 
inner -58" of NGC4826. Contours beyond the dis- 
played field-of-view have been screened as S/N ratio 
is lower than 3 for r>29". Aa and Ad are offsets 
(in arcsec) with respect to the phase tracking center. 
Contours arc -1.75, -0.88, 1.75, 4, 6, 8, 10, 12, 15, 
18, 22 to 34 in steps of 5Jykms -1 beam -1 . The dy- 
namical center is indicated by the star marker. The 
white contour highlights the nuclear region here de- 
fined as CND. l-cr noise = 0.25Jykms -1 beam -1 at the 
center of the image (1-a noise = 0.90Jy km s _1 beam -1 
at the edge of the displayed field-of-view) b) same 
as a) but for the 12 CO(2-l) line. The zoomed view 
shows emission coming from the inner 24". Contours be- 
yond the displayed field-of-view have been screened as 
S/N ratio is lower than 3 for r>12". Contours are - 
2.4 -1.4 2.4 to 19.4 in steps of 2.4Jykms -1 beam -1 . 
Thick contours inside the CND highlight the asym- 
metrical pattern. l-ernoise = 0.36Jykms~ 1 beam -1 (l-cr 
noise = 0.8Jykms -1 beam -1 at the edge of the displayed 
field-of-view). Positions of Arm I and II are indicated. 
Beam-sizes are represented by filled ellipses. 
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Fig. 5. a) Peak brightness 12 CO(1-0) map deprojected 
into the galaxy plane in Log e (r)-$, where r (in arcsec) 
and <& (in degrees) are the polar coordinates. Variable r 
is the galaxy deprojected radius in " relative to the dy- 
namical center and $ is the azimuthal angle in degrees, 
measured counter-clockwise from the western side of the 
major axis. We assume an inclination of z=60° and po- 
sition angle PA=ll2° for the NGC4826 disk (see Sect. 
13.3.11) . The locations of the two m—1 spirals (Arms I and 
II), appearing as two straight ridges, are high-lighted by 
the white-shaded lines, b) Same as a) but for the 12 CO (2- 
1) line. Arm I and the beginning of Arm II are detected. 
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Fig. 6. a) Overlay of the mean- velocity field derived from 
the 1-0 data, in line contours, spanning the range (— 
160 km s -1 , 160 km s _1 ) in steps of 20 km s _1 on the in- 
tegrated intensity map of 12 CO(1-0). Integrated intensity 
as in Fig. 0] starting by 1.75 Jykms -1 beam -1 . Velocities 
are here referred to v sys (LSR)=4:16 km s -1 (thick con- 
tour). Solid black (grey) lines are used for positive (neg- 
ative) velocities. The straight line at PA=112° indicates 
the position of NGC4826's major axis. The star marks 
the AGN position, b) Same as a) but for the 12 CO (2- 
1) line. Integrated intensity contours as in Fig. H start- 
ing by 2.4 Jykms -1 beam -1 . Velocities span the range 
(-120 km s" 1 , 120 km s -1 ) in steps of 20 km s _1 . 
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Fig. 7. a) We show the 12 CO (1-0) position- velocity (p-v) 
diagram along the major axis of NGC 4826 (line contours 
and grey scale range: 0.015, 0.035 to 0.27 Jy beam -1 in 
steps of 0.035 Jy beam -1 ). Velocities have been re-scaled 
to vf^ fl =416 km s -1 and x' offsets are relative to the 
dynamical center, b): Same as a) but for the 12 CO(2-l) 
line in the inner r~15". Contours are the same as in a). 
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Fig. 8. Rotation curve obtained from the combined 1-0 + 
2-1 CO data in NGC 4826 to a radius of R=700 pc. 
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Fig. 9. 12 CO(2-l) position-velocity plot taken along the 
NGC 4826's minor axis. We indicate the position of Arm I, 
and the N/S crossings of the CND lopsided instability. The 
gas kinematics show local signatures of outflow motions at 
these positions. 
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Fig. 10. We display the velocity perturbations of the 
gas flow, projected into the plane of the sky (v per t), 
due to a fast (f2 p =800 km s~ 1 kpc _1 ) (a) and a slow 
(Q p =0 km s _1 kpc _1 ) (b) m=l trailing spiral wave, x and 
y axes are parallel to the major and minor galaxy axes 
(x>0 eastwards, y>0 northwards), and we assume i=-30° 
(i.e., northern side is closer to us). Thin line contours and 
gray scale range from -10, -8 to 16 in steps of 2 km s _1 
(dashed contours for negative values). Isovelocities v=- 
80 km s _1 (dashed line) and v=80 km s _1 (thick line) de- 
fine the orientation of circular rotation in the disk (clock- 
wise). The potential minimum locus is represented by the 
logarithmic spiral. The position of corotation is indicated 
by the dotted circle in the top panel. 
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Fig. 11. We display the overlay of the Paa integrated in- 
tensity image, derived from the HST F187N and F160W 
NICMOS filters (grey logarithmic scale ranges from 
3.1 10 _12 -2.1 10 _11 erg s _1 cm~ 2 pix _1 in the top panel and 
from 7.710 _13 -3.1 10 _11 erg s _1 cm _2 pix _1 in the bottom 
panel; grey scale increases from light to dark), and the 
integrated intensity maps obtained for the 1-0 (top) and 
2-1 (bottom) lines of 12 CO in the disk of NGC 4826 (same 
contours as Fig. EJ) - 
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Fig. 12. The VLA radio-continuum map at 2cm (mostly 
thermal emission) taken from Turner and Ho (1994) (in 
contours from 35% to 95%, in steps of 10% of the max- 
imum=l mjy beam -1 ) is overlaid on the Paa image of 
NGC4826 (in grey scale, as shown). The lopsided star 
formation in the nucleus has a strong maximum offset by 
~16 pc northward from the position of the AGN (marked 
with a square). 
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Fig. 14. Same CO contours as in Fig. II II overlaid on ex- 
tinction Ay derived from the HST broad-band B — I color 
map. Greyscale ranges from Ay=0 to 2. 
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Fig. 15. a(top) Strength of the i-Fourier components 
of the stellar potential of NGC 4826 derived from the 
HST-NICMOS image of the galaxy nucleus. We repre- 
sent the main components-Qi and Q2~and the sum Qt- 
b (middle) The phases of the potential components 
are plotted for i=l,2.c(bottom) Radial variation of the 
average torque-i(i?)-per unit mass. We indicate the radial 
extent of CND, Arm I and Arm II. 



Fig. 16. Estimation of the order of magnitude of the 
Toomre Q parameter (grey scale ranges from 1-3 in the 
top panel and from 0.7-2.5 in the bottom panel) overlaid 
on the integrated intensity maps obtained for the 1-0 (top) 
and 2-1 (bottom) lines of 12 CO in the disk of NGC 4826 
(same contours as Fig. |SJ). 
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Fig. 17. The principal epicyclic frequencies 0, fl-n and 
fl+K in the inner disk of NGC 4826. The pattern speed 
f2 p =1500km s _1 kpc~ 1 has been purposely chosen to lo- 
cate the inner m=l mode of Arm I between its corotation 
and OLR. 



